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Abstract 

Artificial black holes may demonstrate some of the elusive quantum prop- 
erties of the event horizon, in particular Hawking radiation. One promising 
candidate is a sonic hole in a Bose-Einstein condensate. We clarify why Hawk- 
ing radiation emerges from the condensate and how this condensed-matter 
analog reflects some of the intriguing aspects of quantum black holes. 
Keywords: Artificial black holes, Bose-Einstein condensates. 
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1 Introduction 



Picture a Bose-Einstein condensate flowing through a nozzle where the condensate 
exceeds the speed of sound. Suppose that the nozzle is designed such that the 
transsonic flow does not become turbulent. One could build such a nozzle, the 
equivalent of the Laval nozzle [l], out of light, using the dipole force between light 
and atoms to confine the condensate in an appropriate potential. Consider the 
fate of sound waves propagating against the current of the transsonic condensate. 
In the subsonic region sound waves may advance against the flow, whereas in the 
supersonic zone they are simply swept away. No sound can escape the point where 
the flow turns supersonic - the sonic horizon. The transsonic fluid acts as the 
acoustic equivalent of the black hole p], [| . 

An artificial black hole [|J of this kind could be employed to demonstrate some 
elusive quantum properties of the event horizon in the laboratory, in particular 
Hawking radiation H |(| . Hawking predicted that the event horizon emits quanta 
as if the horizon had a temperature given by the gradient of the gravitational po- 
tential. To be more precise, the horizon should spontaneously emit quantum pairs 
where one particle of each pair falls into the hole and the other escapes into space, 
constituting the radiation of the horizon. Both the spectral distribution and the 
quantum state of the emerging radiation are thermal. For solar-mass or larger black 
holes the Hawking temperature is in the order of 10~ 7 K or below, which makes 
the effect next to impossible to observe in astronomy. In the case of sonic holes the 
Hawking temperature is given by the velocity gradient a at the sonic horizon || |3|] , 

*.T=£, (1) 

and the emitted quanta are phonons. A velocity gradient of 10 3 Hz would corre- 
spond to about 1.2 nK temperature. In order to observe such subtle quantum effects 
one should employ the best and coldest superfluids available — Bose-Einstein con- 
densates of dilute gases [f|. 

Some detailed schemes for sonic black holes in Bose-Einstein condensates have 
been investigated theoretically || [K], [lTJ, [12], [13|]. The ultimate design depends on 



experimental details and on the state of the art in manipulating condensates, a 
rapidly evolving field. In this paper we analyze the general aspects of the Hawking 
effect in Bose-Einstein condensates. In the first part of the paper we collect and 
combine the ingredients of the effect, results that are scattered in the literature. 
In the second part we show how the Hawking effect arises naturally within the 



Bogoliubov theory of the elementary excitations in Bose-Einstein condensates WA 



T5|l . F° r the first time, to our knowledge, we connect the quantum physics of the 



event horizon to the behavior of a realistic quantum fluid. 
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2 Sound in fluids 



Consider the propagation of sound in fluids moving with flow speed u. Suppose that 
the flow varies little over the scale of an acoustic wavelength. In this regime we can 
describe sound propagation in geometrical acoustics (the equivalent of geometrical 
optics |L6| or of the semiclassical approximation in quantum mechanics Hl7| ). Sound 
rays follow Hamilton's equations, 

dr du dk du 

~dt = dk ) ~dt = ~dr~' ( ' 

where the dispersion relation between the frequency u and the wave vector k defines 
the effective Hamiltonian u;(r, k). Assume that in each fluid cell 

to' 2 = c 2 k 2 , (3) 

where c denotes the speed of sound and u' refers to the frequency in locally comoving 
frames. In the laboratory frame, uj' is Doppler-shifted, 

u' — uj — u • k . (4) 

In order to see why waves in fluids are related to waves in general relativity, we write 
the dispersion relation (|3]) in a relativistic form. We introduce the space-time wave 
vector 

fc„ = (-w,k) (5) 

and the matrix 

f V = V- 2 ( 1 2 " V (6) 

y \_ u -c 2 l + u <g> u J ( J 

The prefactor Q is an arbitrary non-negative function of the coordinates called the 
conformal factor. In this notation the dispersion relation appears in the relativistic 
form 

g^k„K = 0, (7) 

adopting Einstein's summation convention. Therefore, sound waves experience the 
moving fluid as an effective space-time geometry with the metric g^, the inverse 
matrix of g^, given by 



9p» = ^ 2 



c 2 — u 2 u 
u -1 



The analogy between sound waves in fluids and waves in general relativity turns 
out to be exact for an irrotational fluid |T3| with arbitrary density profile po, flow u 



and speed of sound c, where po, u and c may vary in space and time. The velocity 
potential if and the density perturbations p s of sound obey the linearized equation 
of continuity and the linearized Bernoulli equation |19| 



d tPs + V • (up, + p V^) = 0, (9) 

{d t + u • V)y? + c 2 - = 0. (10) 
Po 
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As a consequence, the velocity potential if of sound obeys the equation || 



d t ^(d t + u-V)<p + V-^[ud t -(c 2 -u 2 )]cp = 0, (11) 



which can be written as the relativistic wave equation 



D v D"<p = -L^v^TU^ = (12) 
V^9 



with the conformal factor Q, in d spatial dimensions, chosen as || 

The assumptions made in order to derive the wave equation (ffl|) are: The fluid [[0| 
is irrotational (1) and isentropic (2). Bose-Einstein condensates naturally satisfy 
condition (1). Condition (2) characterizes the hydrodynamic regime of condensates 
[§]. Here the local pressure depends only on the density and on the temperature of 
the fluid and the quantum pressure is negligible ||. Condition (2) turns out to be 
violated close to the sonic horizon. 



3 Sonic horizon 

Consider the propagation of sound waves in the vicinity of the sonic horizon. Focus 
on the physics in the direction z of the flow at the horizon in a quasi-onedimensional 
model. Assume that the speed of sound in the fluid is constant. The wave equation 
( [LTD reads explicitly 



d t {d t + ud z )ip + d z [ud t - (c 2 - u 2 )d z ] <p = . (14) 
We obtain the general solution 

f dz 

¥ = ¥o( T ± ~t), r± = / —— . (15) 

J c±u 

The t± — t refer to null coordinates in the frame comoving with the fluid |j2U|| . In 
these coordinates sound waves propagate exactly like in homogenous space. In the 
laboratory frame sound waves are accelerated or slowed down by their carrier, the 
moving fluid. An interesting behavior occurs near the horizon, say at z — 0, where 

it = — C + 0£Z . (16) 

The constant a describes the velocity gradient of the condensate at the sonic horizon. 
We see that 

T + = . (17) 

a 

Wave packets localized just before the horizon at z > take an exponentially large 
time to advance against the current. On the other side of the horizon, z < 0, such 
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waves drift equally slowly in the direction of the flow. The horizon at z = marks 
a clear watershed, cutting space into two disconnected regions. In terms of the 
coordinate QT7| ) these regions are characterized by the sign of z^. For stationary 
sound waves with frequency u we get from the general solution (|T5| ) 

(p = Re{<f A z iuj/a e~ lujt } . (18) 

The phase of the wave, (u / a) ln(z / z^,) , diverges logarithmically at the horizon 
where, in turn, the wavenumber k develops a pole, 

k = l(^±) = JL, (19) 
oz \a Zoo J az 

and the wavelength of sound shrinks beyond all scales, 

A = — = \ az. 20 
k 

However, when A reaches the scale of the healing length £ of the condensate || 
(also called the correlation length) the hydrodynamic description of sound in Bose- 
Einstein condensates is no longer valid M. The acoustic theory at the horizon 
predicts its own demise. Similarly, waves near the event horizon of a gravitational 
black hole are compressed beyond all scales. New physics beyond the Planck scale 



may affect the event horizon pi], |22], p3l 



4 Bogoliubov dispersion 

For Bose-Einstein condensates the equivalent of trans-Planckian physics is well- 
known — Bogoliubov's theory of elementary excitations. In the dispersion relation 
(H) we replace the right-hand side by Bogoliubov's famous result [§, |l| [27 1 

,12 „2 7„2 ' - h 



J 2 = c^[l + -J . (21) 
The parameter k c is the acoustic Compton wave number 

with m being the atomic mass, expressed also in terms of the healing length £ || 
(the correlation length). Typically, £ is in the order of 10 _6 m and c reaches a few 
10 _3 m/s in Bose-Einstein condensates (without exploiting Feshbach resonances). 
We calculate the group velocity 

duj , , , 

" = 3k= U + V - (23) 

* = = + (24) 



dk to' \ k 2 c 
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Equation fl23|) shows that the group velocity obeys the Galilean addition theorem of 
velocities. Equation (f^f) expresses the group velocity in the fluid frame, v', in terms 
of the frequency and the wavenumber. The acoustic Compton wavenumber, k c , sets 
the scale beyond which v' deviates significantly from c. For large wavelengths, sound 
is communicated by atomic collision, and the product of the condensate's density 
and the atomic collision strength gives mc 2 ||. For wavelengths comparable with 
or shorter than the healing length, the interaction-free Schrodinger dynamics of the 
atoms dominates the transport of excitations. Perturbations of the free wavefunction 
travel with infinite velocity. So the acoustic Compton wavenumber, k c , characterizes 
the crossover between the speed of sound and the infinite speed of perturbations of 
free matter waves. 

Close to the sonic horizon, the wavenumber ([19]) increases dramatically, and, 
in turn, the effective speed of sound v' grows. The horizon, defined as the place 
where the fluid exceeds the speed of sound, seems to dissolve like a mirage. Nature 
appears to prevent the existence of an event horizon. However, we show in Sec. 6 
that the horizon still exists, but at a less well-defined location and for a particular 
class of elementary excitations only. As long as k 2 is much smaller than k 2 we get 
the acoustic relation 

k — . (25) 

u±c 



For the other extreme, where k is much larger than k 2 , one finds [28, 29 



k ~ ±2k c ^uyc 2 - 1 + . (26) 



c 2 — u 2 



Consider the turning points zq, the points where the group velocity of sound (|23|) 
vanishes. If the acoustic dispersion relation (|3]) were universally valid the horizon 
would be the turning point. Therefore \z\ < \z \ does indicate the spatial scale of the 
trans-acoustical range around the horizon, which defines the spatial derealization 
of the horizon. To proceed we recall that elementary excitations are small pertur- 
bations of the condensate. Their energies Hu ought to be much smaller than the 
mean-field energy of the condensate, which is in the order of mc 2 (with c being the 
speed of sound). Therefore, 

(27) 



mc z 



We expand the solution z of v = as a power series in e 1 ^ 3 and find, to leading 
order, three turning points in the complex plane given by 



C3 3/ /£\2/3 



a2 \2J 

Far away from the horizon we may characterize the four fundamental solutions of the 



dispersion relation (21) combined with the Doppler shift (|J) by their asymptotics 
( |25"| ) or (]26f) . However, close to a turning point geometrical acoustics alone does 
not provide a good description of wave propagation anymore. The turning points 
may cause scattering. The connections between elementary excitations across the 
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horizon must be examined with care. Extending z to the complex plane represents 
an elegant way of analyzing this connection. We find in the Appendix that the 
acoustic relation (|T^) remains valid on either the upper or the lower half of the 
complex plane. 



5 Bogoliubov modes 

Elementary excitations are perturbations of the condensate, ripples on the macro- 
scopic wave function ipo of the condensed atoms. The excitations constitute the 
non-condensed part of the atomic gas. To describe elementary excitations, the total 
many-body field operator if) of the atoms is split into two components, the con- 
densate with the mean- field wave function tp and the non-condensed part. The 
mean-field wave function comprises the density profile p an d the flow u, as 

^o = ^oe iS \ u=-VV (29) 

m 

The atomic field operator ip is split into the condensate and the non-condensed part 
according to the relation 

$ = + e iSo 4> . (30) 
The non-condensed part consists of Bogoliubov modes u v and v v J8], [b|, [L5| , 



+ v t°l) ■ ( 31 ) 



The u v and v v are subject to the Bogoliubov-deGennes equations M, 14, 15|. If one 



requires that the Bogoliubov modes satisfy the orthonormality relations 

{u* u u u > - v* v v v >) dz = 5 w i , (32) 

{UyVyl — VyUyl) CiZ = 0, (33) 



then the d u and aj, obey the commutation relations of Bose annihilation and creation 
operators, as a consequence of the fundamental commutator of atoms with Bose 
statistics 

i>(z),ft{z')\ =5{z-z'). (34) 

Each pair of u v and v v characterizes the spatial shape and the evolution of an 
excitation wave, and the Fock space of the a v and aj, spans the state space of the 
excitation quasiparticles. 

To see how the Bogoliubov modes are related to sound waves, we write down the 
macroscopic wave function of the condensate combined with one of the excitations, 

i) = ip + e tSo (u u + vl). (35) 

We represent ip as 

^ = A /pe i5 , p = Po + Ps , S = S + s, (36) 
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where p s denotes the local density of the sound wave and s is proportional to the 
velocity potential 

ip=—s. (37) 
m 

Assuming that p s and s are small perturbations we get 

u v + v* u = ^ (jjL + is^J . (38) 
Assuming further that u v and v v are stationary waves with frequency uj, we obtain 



from our solution (18) of the hydrodynamic sound-wave equation and from the 



linearized Bernoulli equation (|TUp the Bogoliubov modes 

v v ~ AJ^-^z^e-^. (39) 
\2az n / 

These asymptotic expressions are valid as long as the elementary excitations are 
sound waves with wave number (0). We show in the Appendix that this is the case 
sufficiently far away from the turning points and on either the upper or the lower 
half plane. Here the term co/(2az) in the expressions (|39|) is always small compared 
with mc/H. For real and positive z we have 

( — £f w l a = e -(±27ruV°0 z iu/a _ 

The ± sign refers to the two ways in which we may circumvent the trans-acoustic 
region, on the upper (+) or on the lower (— ) half plane. Modes with the acoustic 
asymptotics (|3^) throughout the upper half plane are suppressed on the left side of 
the horizon and for positive frequencies uj and enhanced for negative to. Modes with 



the asymptotics (p9|) on the lower half plane show the opposite behavior. 



6 Negative energy 



Bogoliubov modes are normalized according to the scalar products fl32|) and (|33|). 
Let us calculate the norm of the modes (|39|) of the sonic hole. The energy parameter 
( p?| ) is small and so is the extension of the trans-acoustic region around the horizon, 
measured roughly by the location of the turning points ( p8|) . Consequently, we 
can neglect the trans-acoustic contribution to the normalization integral (p2). We 
approximate the Bogoliubov modes by their asymptotic expressions (j3~9|), utilize the 
relation ([40]), and get 



{u* v u v , - v* v v v .) dz ~ \A U \ 2 ( f + [ + °°)^ z ^'-y-dz 

\J~oo Jo ) az 

= \A U \ 2 (1 - e -(±W°0) ^cukJiu' - u) . (41) 

If we choose modes with the asymptotics ( p^) on the upper half plane the norm is 
positive and the u v and v v may serve as proper Bogoliubov modes. We find the 
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normalized amplitude 



A v = [(l - e" 2 ™/") 4nuk c ] 1/2 . (42) 

Remarkably, the Bogoliubov norm is positive also for modes with negative frequen- 
cies. The Bogoliubov-deGennes equations have a well-known symmetry |14|, |T5| : 
If the (u u ,v v ) are solutions then the complex-conjugated and interchanged modes, 
(f*,u*), are solutions as well. Yet the norm of the conjugate modes is the nega- 
tive norm of the original (u u ,v u ). In contrast, sonic black holes generate negative- 
frequency modes with positive norm, which is the unusual feature that gives rise 
to the acoustic analog of Hawking radiation ||. We also see how the mentioned 
symmetry of the Bogoliubov-deGennes equations [14 , H| appears in our case. If we 



chose the u u and v v with the asymptotics fl39|) on the lower half plane we would get 
negative normalization integrals ( 4l|) for both positive and negative frequencies. 



The negative frequencies of the positive-norm Bogoliubov modes give rise to 
negative energies in the Hamiltonian of the elementary excitations, 

H = J Huj (a+a+ - aLa-^j N du . (43) 

Here and later the subscripts ± refer to positive and negative frequencies, respec- 
tively, and N denotes the density of modes. We see that there is no natural ground 
state of the elementary excitations. In practice, of course, the spectrum is limited by 
the requirement ([27|) that the energies of the excitations ought to be much smaller in 
magnitude than the condensate's mean-field energy. We note that the Hamiltonian 
(|43"1) is invariant under the Bogoliubov transformations 

a' ± = a± cosh £ — sinh £ , 
u'± = u± cosh £ — d* sinh £ , 

v' ± = v± cosh £ — u* sinh £ . (44) 

with an arbitrary real parameter £. In the case we choose 

tanh£ = e-™ /a (45) 

we get a new set of modes, Eq. ([39]), with 

A' ± ~ 9 {±z) (Aixujk c )- 1/2 . (46) 

The function indicates that the primed modes appear on either the left or on 
the right side of the trans-acoustic region. Therefore, despite the trans-Planckian 
problem, a sonic horizon exists, but at a less well-defined location, within \z\ < \zq\, 
and the horizon applies to a particular set of modes only. 
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7 Hawking effect 



The Bogoliubov transformations (fHj) relate one set of quasiparticles to another one, 
both representing perfectly valid energy eigenvalues, yet their quasiparticle vacua 
differ (the states |0) or |0') that are annihilated by a± and a' ± , respectively). This 
ambiguity has direct physical consequences, because the cloud of non-condensed 
atoms [26| depends on the vacuum state of the elementary excitations, 



(0\ft4>\0) = J (K| 2 + \v„\ 2 ^Ndu, (47) 
(O'|0t0|o') = / (\v' + \ 2 + \v'_\ 2 ^Nduj 



± <0| 0+0 |0> . (48) 

In general relativity the notion of the vacuum is observer-dependent. For example, 
the vacuum of empty space in Minkowski coordinates appears as a thermal field to 
accelerated observers |5], ||, [H]]. In the case of the black hole, the gravitational col- 
lapse has created a state of quantum fields that an inward-falling observer perceives 
as vacuum, yet an external observer sees as thermal radiation, Hawking radiation 
[|> S 01- I n our case, the equivalent of the gravitational collapse, the formation of 
the sonic horizon, chooses the quasiparticle vacuum, if the trans-sonic velocity pro- 
file has initially been created from a condensate without a horizon. Such a process 
must be sufficiently smooth to keep the condensate intact. 

To analyze the quasiparticle vacuum, we use the Heisenberg picture of quan- 
tum mechanics where observables evolve while the quantum state is invariant. We 
describe the initial (and final) vacuum state with respect to one set of continuous 
modes given before the formation of the horizon. In the Heisenberg picture these 
modes evolve. We sort the initial modes into left- and right-moving modes that, 
close to complex oo, are analytic on the upper or on the lower half plane, respec- 
tively, because here exp(ikz) converges for positive k on the upper and for negative 
k on the lower half plane. The upstream modes we are interested in stem from 
right-moving modes. The formation of the sonic horizon, a smooth process, cannot 
fundamentally alter the analyticity of the vacuum modes. In particular, the process 
can never create non-analytic modes of the type expressed in Eq. (|46|) . Consequently, 
the initial quasiparticle vacuum assumes the analytic modes of Eqs. ( PU|) and ([K|). 

Given this vacuum state, we determine the quantum depletion of the condensate. 
We write the density of the non-condensed atoms in terms of the primed Bogoliubov 
modes. As we have seen, these modes describe the set of elementary excitations that 
exhibit the sonic horizon. We find that 

(O|0t0|O) = p(|2|), p{z)= /[(|<| 2 + K| 2 )nH + K| 2 liV^. (49) 



Here n(u) denotes the average number of non-condensed atoms per excitation mode, 

, . 1 1 
niuj) = — — ; = . ,, = . (50) 

v 1 e 2-KLu/a _ ]_ e huj/k B T _ j v 1 
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The non-condensed atoms are Planck-distributed with the Hawking temperature 
(0). Therefore, as soon as the condensate flows through the nozzle, breaking the 
speed of sound, a thermal cloud of atoms is formed. This effect is the signature 
of Hawking radiation for sonic holes in Bose-Einstein condensates. The thermal 
cloud due to the Hawking effect should be observable when the initial temperature 
of the atoms is below the Hawking temperature. On the other hand, one could 
also regard the Hawking effect in the condensate as the quantum depletion ([|7]) of 
atoms at zero temperature with respect to the analytic modes ([42]) that transcend 
the horizon. This feature reflects the ambiguity of the vacuum in general relativity. 
Using techniques for measuring the population of Bogoliubov modes |32], ^3|, one 
could perhaps demonstrate the ambiguity of the vacuum in the laboratory. 

Finally we note that within our model sonic black holes are stable, provided 
of course that the transsonic flow is not plagued by hydrodynamic instabilities. 
In reality, elementary excitations interact with each other, giving rise to what is 



known as Landau-Beliaev damping [34, 35, 36]. Since a sonic horizon does not 



have a ground state, this damping mechanism will lead to the gradual evaporation 



of the condensate. Therefore, Landau-Beliaev damping ||34| , pq , |36fl plays the role 
of black-hole evaporation. It is tempting to turn matters around and to approach 
cosmological problems from the perspective of condensed-matter physics |37|, |38|, |3T 
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Appendix 



In this Appendix we examine the asymptotics of the Bogoliubov modes on the 
complex plane. We use the semiclassical approximation pi HO] 



u u = U u exp yi J kdz — iuit 

v u = K,exp ( i / kdz — iujt I . (51) 



The wavenumber k should obey Bogoliubov's dispersion relation (21) including the 



Doppler shift (|). We obtain four fundamental solutions of this fourth-order equa- 
tion. Figure 2 shows the three branches of k that are relevant in our analysis. The 
amplitudes U v and V v obey the relation |30], f4(J 

d z (U 2 U -V u 2 )v = 0, (52) 

where v denotes the group velocity (Effi) of the elementary excitation. Equation (|52|) 
formulates the conservation law of the quasiparticle flux for stationary states [[4pJ 
if z is real, where U 2 — V 2 gives \U U \ 2 — \V U \ 2 up to a constant phase factor. The 
relation (B^) can be extended to the complex z plane and to complex frequencies u 
PHI . Equation (|52|) implies that the amplitudes U u and V v diverge close to a turning 
point z where v vanishes. Consequently, the semiclassical approximation breaks 
down at the turning point [|I7|. The turning point causes significant scattering, i.e. 
the conversion of one mode with a given k into two modes, one with wavenumber k 
and the other one with a different wavenumber that satisfies the dispersion relation 
as well. At the turning point the two branches coincide. To prove this, we regard 
for a moment z as a function of k at constant frequencies uj = ujq, defined implicitly 
by Eqs. (§, flU) and (0). We get 

v = - = u + m { Uo -vk) = -ak-. (53) 

We see that the function z(k) reaches extrema at the point Zq where v vanishes, i.e. 
at the turning point. Close to the z we find after some algebra by expanding 
z into a power series in e 1//3 , 



3{k-k ) 2 , 3r — 
z-z , k = V^Ai , 54 

to leading order. Consequently, each turning point connects two branches of the 
wavenumber k. In general the mode conversion occurs near specific lines in the 
complex plane, called Stokes lines in the mathematical literature |41| . Stokes lines, 



originating from the turning point zq, are defined as the lines where the differences 
between the phases J kdz of the two connected k branches is purely imaginary. Here 
one of the waves is exponentially small compared with the other. We obtain from 
Eq. fl5~4]) that the difference between the two branches is proportional to the square 
root of z — zq. Consequently, the phase difference is proportional to (z — z ) 3 ^ 2 , 
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giving rise to three Stokes lines from each turning point z , as in the traditional case 



of Schrodinger waves in one dimension [17, 42|. Figure 3 shows the Stokes lines for 



the three turning points close to the horizon and for the branch cuts of k chosen in 
Fig. 2. 

At a Stokes line the phase difference between the two connected branches is 
purely imaginary. One of the waves exponentially exceeds the other and, within 
the semiclassical approximation, the smaller wave is totally overshadowed by the 
larger one, if the larger wave is present. In general, the Bogoliubov modes consist 
of a superposition of the four fundamental solutions that correspond to the four 
branches of the dispersion relation 

u v c A u A -\- c B u B -\- c c ii c -\- C D U D , 

v v = c A v A + c B v B + c c v c + c D v D . (55) 



The u A and v A refer to the k branch where k obeys the asymptotics (19), i.e. where 
the wavenumber satisfies the dispersion relation ([$]) of sounds in moving fluids, tak- 
ing into account the Doppler detuning (f|) and where k corresponds to an upstream 
wave. We call such Bogoliubov modes acoustic modes. When crossing a Stokes line, 
the exponentially suppressed solution may gain an additional component that is pro- 



portional to the coefficient of the exponentially enhanced solution |42|. If we wish to 
construct Bogoliubov modes where only the exponentially smaller component exits 
in the vicinity of a Stokes line we must put the coefficient of the larger one to zero. 
In Fig. 3 the pairs of letters indicate which branches are connected by the lines, and 
the first letter identifies the exponentially dominant branch. The picture shows that 
with the choice of branch cuts made we can construct a Bogoliubov mode that is 
acoustic on the upper half plane. Trans-acoustic physics is confined to the lower half 
plane. On the other hand, if we chose other branch cuts of k we may get Bogoliubov 
modes that are acoustic on the lower half plane and trans-acoustic on the upper one. 
Therefore, according to Eq. fl40|) , the choice of the k branch determines whether a 
Bogoliubov mode is larger or smaller beyond the sonic horizon at the real axis, for 
z < 0. Branch cuts of k are fairly arbitrary. Given the Bogoliubov mode on the right 
side of the horizon, we cannot predict within the semiclassical approximation the 
amplitude of the mode on the left side. Therefore, the two sides are causally discon- 
nected. Within the semiclassical approximation, the horizon is a genuine horizon, 



despite the acoustic analog of the notorious trans-Planckian problem [21, 22 
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Figure 1: Schematic diagram of a sonic horizon. A fluid is forced to move through a 
constriction where the flow speed u becomes supersonic (dashed line). The constric- 
tion may be formed by the walls of a tube or, if the fluid is an alkali Bose-Einstein 
condensate, by a suitable trapping potential. The picture shows a Laval nozzle 
where a supersonic fluid is hydrodynamically stable. 
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Figure 2: Wavenumbers k of elementary excitations around a sonic black- hole hori- 
zon, analytically continued onto the complex plane. The figure shows three roots 
of the dispersion relation [h 2 k 2 /(2m) + mc 2 ] 2 — h 2 [uj — (— c + az)k] 2 = m 2 c 4 for 
ou — 0.1 (mc 2 /h) and a = 0.5 (mc 2 /h), illustrating the branch cuts of k. The top 
row displays the wavenumber of a sound wave that propagates against the current. 
The picture indicates the characteristic u/(az) asymptotics away from the branch 
points. The two lower rows display two trans-acoustic branches of k. The fourth 
root of the dispersion relation is not shown, because it corresponds to the trivial 
case of sound waves that propagate with the flow. 
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Figure 3: Stokes lines of elementary excitations at a sonic black hole (dotted lines), 
given the choice of branch cuts made in Fig. 2. The pairs of letters indicate which 
branches of the superposition (52) are connected by the lines. The first letter of each 
pair identifies the exponentially dominant branch. We construct a Bogoliubov mode 
that is acoustic (component A) on the upper half plane by putting the coefficient 
c c to zero on the C,A Stokes line originating from the left turning point. 
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